PROGNOZ-3S

1. PROGNOZ-3S

Решается задача по трехмерной расчетной схеме. Главная (управляющая) программа MAIN составляется для каждого варианта прогноза. В программе MAIN вводятся исходные массивы начальной температуры и некоторые другие исходные данные. Размерности массива, индивидуальные для каждого варианта прогноза, задаются в операторах описания.

Из программы MAIN происходит обращение к подпрограмме WNS. В ней определяются коэффициенты А, В, и С для вычисления количества незамерзшей воды в каждом слое грунта, по формуле где U£0 - температура грунта. Для контроля выдается на печать таблица получаемых значений W_w от U при.

Из подпрограммы WNS получены массивы А(N), В(N), С(N) передаются в MAIN.

В подпрограмме GRUND определяются одно-, двух- или трехмерные температурные поля в заданном массиве грунтов. Здесь используются подпрограмма-функция HF и подпрограмма FIND. HF (N,U1) вычисляет энтальпию 1 м³ грунта в N-ом слое при температуре U1<URK (критическая температура фазовых переходов). Для расчета функция HF использует физико-механические и теплофизические свойства грунтов, введенные подпрограммой BLOCK DATA.

Подпрограмма FIND (HN1, N1,UL1) определяет температуру UL1 в слое грунта с номером N1 по энтальпии HN1, если HN1<H1(N1), т.е. еслисоответствует энтальпии мерзлого грунта. FIND использует общую область TABL, в которой содержатся таблицы значений энтальпии по температуре от 0 до -50°С по всем грунтовым слоям, вычисленным с помощью функции HF.

Подпрограмма ISTO (i, j, K, T, F) учитывает влияние источников тепла и охлаждающих установок, находящихся внутри массива грунта. Результат работы подпрограммы присваивается простой переменной F, где F - тепло, выделяемое внутренними источниками (стоками) в единицу времени (час) в фиксированном элементе (i,j,K) грунтового массива в момент времени T.

Алгоритм подпрограммы может быть составлен как по экспериментальным данным, так и по теоретическим расчетам применительно к конкретному случаю.

Для обращения к подпрограмме ISTO в подпрограмме GRUND (рекомендуемое приложение 2) следует заменить оператор А=0.0 с номером 0000399 на обращение к подпрограмме ISTO:

Источник: РСН 67-87: Инженерные изыскания для строительства. Составление прогноза изменений температурного режима вечномерзлых грунтов численными методами